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Motion of two gravitating spherical stellar shells around a massive central body is considered. 
Each shell consists of point particles with the same specific angular momenta and energies. In the 
case when one can neglect the influence of gravitation of one ("light") shell onto another ("heavy") 
shell ("restricted problem") the structure of the phase space is described. The scaling laws for the 
measure of the domain of chaotic motion and for the minimal energy of the light shell sufficient for 
its escape to infinity are obtained. 



PACS numbers: 95.10.Fh; 98.10.+Z 

Chaotic motion of a system can result in 
crossing a border in the phase space sepa- 
rating the area where disintegration of the 
system is possible. In the present paper this 
phenomenon is studied in a system of two 
gravitating stellar shells oscillating around 
a massive central body. Similar phenomena 
take place in planetary and molecular dy- 
namics. 



I. INTRODUCTION 

Dynamical processes around supermassive black 
holes (SBH) in quasars, blazars and active galac- 
tic nuclei are characterized by violent phenomena, 
leading to formation of jets and other outbursts. 
Here we consider a model problem for a shell out- 
burst from a SBH surrounded by a dense massive 
stellar cluster due to a pure gravitational interac- 
tion of shells oscillating around SBH. 

Investigation of spherical stellar clusters using 
shell approximation was started by Henon |^. In 
this approximation a stellar cluster is considered 
as a collection of spherical stellar shells. Each shell 
consists of stars with the same specific angular mo- 
menta and energies. In the process of motion the 
angular momenta are conserved while energies are 
changing, yet taking the same values for all the 



* Electronic address: bamiv@sai.msu.ru| 
tElectronic address: gkogan@mx.iki. rssi.ru| 
t Electronic address: aneishta@m x.iTti.rssi.rul 
''Electronic address: i.belinski@icra.itl 



stars in one shell. The shell itself stay spherical. 
Shell approximation has been successfully applied 
for investigation of stability , violent relaxation 
and collapse , leading to formation of a sta- 

tionary cluster. Investigation of the evolution of a 
spherical stellar cluster taking into account various 
physical processes was done on the base of the shell 
model in the classical series of papers of L. Spitzer 
and his coauthors IIHHIllillilEllIalll, 
see also [T^ . 

Numerical calculations of a collapse of stellar 
clusters in theshell approximation 0, have 
shown, that even if all shells are initially gravita- 
tionally bound, after a number of passages through 
each other some shells obtain sufficient energy to 
be thrown to infinity. 

Here we consider a simplified problem of motion 
of two massive spherical shells, each consisting of 
stars with the same specific angular momentum 
and energy, around SBH. 

Development of chaos during the motion of two 
gravitating intersecting shells was found first by 
Miller and Youngkins [IjI in an oversimplified case 
of purely radial motion and reflecting inner bound- 
ary. In 0, chaotic behavior was found in a more 
realistic model considering stars with the same 
specific angular momentum and energy dispersed 
isotropically over the spherical shell, and each star 
moving along its ballistic trajectory in the gravi- 
tational field of the central body and shells. In the 
current paper we consider the same model mainly 
in the case when one can neglect the influence 
of one ("light") shell onto the motion of another 
("heavy") shell. In this simplified case we describe 
a structure of the phase space and obtain the scal- 
ing laws for the measure of the domain of chaotic 
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motion and for the minimal energy of the hght shell 
sufficient for its escape to infinity. 

We follow the approach which was used in |20| 
for analysis of chaotic dynamics and escape to in- 
finity in the planar restricted circular three body 
problem. This approach consists of constructing 
asymptotic formulas for Poincare return map and 
identification of regions of regular and chaotic dy- 
namics in the phase space of this map. Region of 
regular dynamics is the region where Kolmogorov- 
Arnold-Moser (KAM) theory (see, e.g. [1|3|) guar- 
antees existence of many invariant curves of the 
map. These curves are barriers preventing the es- 
cape to infinity. Region of chaotic motion is the 
region where the phase expansion criterion (see, 
e.g. |33|) is satisfied. In this region chaotic diffu- 
sion allows the " light" shell to gain positive energy 
and then escape to infinity. One special property of 
the system under consideration is that for an open 
set of the problem's parameters, in the phase space 
there is a region where the Poincare return map is 
continuous but not smooth (it has a singularity of 
the square root kind). In this region neither KAM- 
theory nor phase expansion criterion are applica- 
ble. Similar situation, though with a map that is 
discontinuous in a certain region in the phase space 
occur in the problem of motion of three gravitat- 
ing parallel impenetrable sheets [l^- For general 
analysis of maps with discontinuities see Q and 
references there. 

The dynamics of a system of many gravitating 
spherical shells was studied from the viewpoint of 
statistical mechanics in |0, 0, where, how- 
ever, the escape to infinity was not studied because 
a bounded system was considered. An extensive 
study of escape to infinity in the rectilinear three 
body problem was performed in . 



II. 



A RESTRICTED PROBLEM. 



Consider the motion of two stellar shells with 
masses mi and m2 around the central body of mass 
M. Let Ji and J2 be specific angular momenta 
of shells. Consider the case when 7712 <C mi, so 
that one can neglect the infiuence of the second 
shell onto the motion of the first shell. We call 
this problem a restricted problem by the analogy 
with the restricted three body problem in celestial 
mechanics (the latter problem describes the motion 
of an asteroid under the action of Sun and Jupiter) . 

In what follows the first shell will be called the 
heavy one, and the second shell will be called the 
light one. 



In the restricted problem, motion of the heavy 
shell is described by Hamiltonian function 
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G{M + mi/2) J\ 



ri 



Here pi is the momentum canonically conjugated 
to the shell radius ri. Therefore, the change of ri 
in time is the same as in the Kepler problem for 
a satellite with mass 1 and angular momentum Ji 
and a central body of mass {M ~\-mi/2). For Hi = 
hi < Q the behavior of ri is described by usual 
formulas of the elliptic Keplerian motion: ri = 
ri(/i, /ii, Ji); here li is the mean anomaly of the 
satellite, i.e. the angular variable for which 

'1 = ^' Ml = G(Af + mi/2), ai^--^. 
ai ^"1 

(2) 

Motion of the light shell is described by Hamilto- 
nian 



H2 
V2 



^2 IT/ 



if r2 <ri, 



0, ' Z ^ ' 1, /ON 

Vb — Gmi/r2, if r2 > ri. 



where 



GM 
r2 



Jl 
2t\ 



Here p2 is the momentum canonically conjugated 
to radius T2 of the light shell. In the Hamiltonian 
li2 one should consider ri as a given function of 
time. Therefore, in the restricted problem motion 
of the light shell is described by a Hamiltonian sys- 
tem with one and a half degrees of freedom (one 
degree of freedom plus explicit dependence of the 
Hamiltonian on time). This model applies just 
as well to the case where the light massless shell is 
replaced by a massless point. 



III. POINCARE RETURN MAP FOR 
THE RESTRICTED PROBLEM. 

Consider the motion of the light shell. At the 
time moments t, when r2 has a local minima 
on the trajectory, we mark on the cylinder ^ = 
{{h,(j)) : —00 < h < 00, <j) mod 27r} the values of 
the light shell Hamiltonian H2 and of the heavy 
shell mean anomaly li : {H2,li) = {h,(j)). Thus 
we obtain on $ a sequence of points . . . , (/i^, (/)*), 

...,(/^-^r'),(/l^^o),(/^^^l),...,(/l^</.^),.... 
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This sequence terminates at the right if the sheU 
escapes to infinity. It can be terminated at the 
left as well, but the terms of this sequence with 
numbers smaller than certain negative number do 
not have physical meaning. The points of the con- 
structed sequence are mapped to each other under 
the action of Poincare return map defined as fol- 
lows. Consider the motion of the light shell that 
starts at time moment t = t with P2 — 0. Denote 
initial values of H2 and h as h and (p. The ini- 
tial value of r2 satisfies the equation H2 = h. Let 
this initial value r2 be the smallest of two roots 
of equation H2 = h. Therefore, at the time mo- 
ment t the value of r2 on the trajectory has a 
local minimum. Denote t' as the first time mo- 
ment after t such that at this time moment again 
P2 = and r2 on the trajectory has a minimum, 
provided that such a time moment exists. Denote 
as h! and (j)' values of H2 and li corresponding 
to t — t' . By definition, Poincare return map 11 
is given by the formula n(/i, (/>) = {h',(f)'). Then 
(/j(fe+i)^0(fe+i)) ^ n(/i('=), The map H is 

defined on the part of the cylinder $ where values 
of {h, 4>) allow the shell to make at least one com- 
plete oscillation. The map 11 preserves the area 
dhdcj) on $; this follows from preservation of the 
phase volume by a Hamiltonian system (see, e.g. 
Refs. 

Together with map 11 we will use modified re- 
turn map n defined as follows. Let t,t',(j),<j>' 
be the same as before. Denote by h the value 
of H2 at the last preceding to t time moment 
when r2 has a maximum. Denote by h' the 
value of H2 at the first time moment after t 
when r2 has a maximum. Define map 11 by 

the formula fl(ji,(l?j — The map 

n is defined on the part of the cylinder $ — 
I (ji, (/)^ : —00 < h < 00, (j) mod 27r| where h < 

0. The map IT preserves the area {^dh/dfij dhdcf) 

on $. For a given {h},(lf^ the map ft gener- 
ates on $ the sequence of points . . . , (/i'*, 0'*), . . . , 

(/i-i, r (^°, 4'°)Xh\ ct^^). . . . , 0"), . . . such 

that (/i('=+i',</)(''-+i)) = ft(ftW,0Wy At smah 

values of h the time interval between t and t' is de- 
termined in the main approximation by the value 

h': t' -1=21^112 I {^2h!^''^ , ^i2=G{M + mi). 

This is the reason why for mall h map 11 is more 
convenient for the analysis than 11. In what follows 
we omit for brevity the symbol """ over h, 



In what follows we consider the case when the 
heavy shell's mass is much smaller than the mass 
of the central body: mi = eM, < e <C 1. In this 
case map ft has the form 

h' ^h + efih,(j),e), (4) 



^' = + 27r^ (^^^ +e9{h',^',e). (5) 

Functions /, g can be extended as continuous func- 
tions up to the argument value h ^ 0. The value 
ef {h, (j>, e) is the sum of jumps of the light shell's 
potential energy due to passage of shells through 
each other between two consecutive time moments 
when the light shell radius takes maximal values. 
To find of e/ {h, cj), 0), one should calculate the sum 
of the energy jumps under the assumption that mo- 
tion of the light shell is not affected by the heavy 
shell (i.e. this motion is Keplerian one with total 
energy h). 

Depending on the properties of the shell motion 
at £ = 0, one can define three main types of h 
values. 



1. Type A. 

At all initial phases (p of the heavy shell mo- 
tion, the light shell on the period of its motion 
passes through the heavy shell at non-zero rela- 
tive velocity. (If h — 0, i.e. motion of the shell 
stars is parabolic, then in this and following defi- 
nitions one should consider as a period of the light 
shell motion the whole infinite time interval be- 
tween arrival of the shell from infinity and its de- 
parture to infinity.) In this case during a period 
of the light shell motion the number of passages 
of shells through each other is the same for any 
initial phase of the heavy shell motion. Therefore, 
the number of terms in the sum for / (/i, (j), 0) is 
the same for any (p. At small enough e, the sum 
for / {h, (p, e) has the same number of the terms. 
Therefore, function f{h,(p,e) is analytic in {h,(p), 
and it can be continued analytically in some neigh- 
borhood of a circle around h ^ 0. 



2. Type B. 

There exists certain initial phase 0, of the heavy 
shell motion such that at some time moment the 
shells superimpose with zero relative velocity and 
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non-zero relative acceleration, and absolute veloc- 
ities of both shells at this moment are different 
from zero. At this moment trajectories of the 
shells have quadratic tangency. When (j> deviates 
from 0* in one direction, the superimposing of the 
shells disappears. When (p deviates from 0* in 
another direction, this superimposing decomposes 
into two close shell crossings (say, light shell sur- 
passes the heavy one and then almost immediately 
the heavy shell surpasses the light one). In this 
latter case two new terms appear in the sum for 
f{h,(j),0); these terms are close by absolute val- 
ues, and they have differ ent signs . The sum of 
these terms is of order ^/{(j) — (/)*|. Therefore, in 
this case function / {h, (j), 0) has at = the sin- 
gularity of type ^/(j) — 4>^, or y/cj)^: — (j). At small 
variations of h this singularity persists. Func- 
tion / (/i, (f>, e) has the same singularity at small 
e. Value r2 that corresponds to time moment of 
tangency of shell trajectories is uniquely defined. 
Indeed, at the moment of tangency there should 
be ri = r2, hi + Jll2r\ = h + J|/2r| (because 
velocities of the shells are equal to each other) 
and Jf/2rf ^ J| /2r2. These conditions define a 
unique value of r2 ■ 

3. Type C. 

For any initial phase (j) of the heavy shell mo- 
tion the shells never superimpose. Then at small 
enough e motion of the light shell is also Keplerian. 

There is also exceptional intermediate type of 
h for which the shells can meet with zero veloci- 
ties (i.e. one of two extremal values of radius of 
one shell coincide with an extremal value of radius 
of another shell), but with nonzero relative accel- 
eration. Another exceptional case is h = hi and 
Ji = J2. In this case for some initial phase (j) of 
the heavy shell the radii of two shells coincide all 
the time. 



IV. DYNAMICS OF SHELLS IN THE 
RESTRICTED PROBLEM. 

Let us consider first the case when the value 
/i = is of type A. Consider in cylinder $ the 
annulus a < h < such that all values h on this 
annulus are of the type A (here a is a constant, 
a < does not depend on e). In this annulus the 
map n is analogous to " Kepler map" , which was 
introduced in (2^ . Let us repeat for map II the 



analysis from (see also 0| ). The trajectories 
of several initial points under the action of map 
n for this case are shown in Fig. ^ The proper- 
ties of the map to large extent are determined by 
the phase expansion coefhcient {dcf)' / d(f)j^. In our 
case relations (0) imply that 

In the region where {d(f>'/d<j>) ^ 1, i.e. h < 
—KihiE^/^ with a big enough positive constant 
i^i, map n satisfies the assumptions of KAM- 
theory. According to KAM theory, this region is 
filled up to a residual set of small measure by in- 
variant curves of map 11 that enclose cylinder $ 
(i.e. values of </> along any of these curves cover 
all segment [0, 2tt] mod 2tt) and are close to curves 
h = const. The residual set of small measure is 
filled up to another residual set of much smaller 
measure by invariant curves that do not enclose 
cylinder $ (i.e. values of along each such curve 
cover only a part of segment [0, 27r] mod 2tt). 
These invariant curves are organized into " islands" 
corresponding to different resonances. For phase 
points on invariant curves the motion of the light 
shell is quasi-periodic (in particular case, periodic). 
For a phase point not on invariant curves the mo- 
tion of the light shell can be chaotic. However, 
because the trajectory of the phase point on re- 
mains locked between invariant curves, the light 
shell can not acquire the positive energy and es- 
cape to infinity. The invariant curves with largest 
h value that can be constructed by KAM theory 
have h ~ — e^/^, the variation of h along any such 
curve is also of the order of e^/^. Fig. [^represents 
results of numerical experiment in which for several 
values of e we looked for the invariant curve with 
the maximal value of h enclosing cylinder $. The 
least mean square fit gives the scaling h ~ — e° '*^. 

In the region where {d(j>' /d(j>) ^ 1 small changes 
in (f) lead to large changes in 0'. Therefor, values 
of and (j)' can be considered as statistically inde- 
pendent according to to phase expansion criterion. 
This is the region where -K^^hie"^^^ < h < Q, 
with K2 being a large enough positive constant. 
In this region, variations of the light shell energy h 
during one iteration of map 11 can be treated as in- 
dependent random values. As a result the changes 
of the energy along the phase point trajectory are 
of diffusional type. Apparently, as a result of this 
diffusion almost all phase points in the region un- 
der consideration acquire positive energy, i.e. light 
shell escape to infinity. Numerical experiment (see 
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Figure 1: Trajectories of five initial points (3 chaotic, 2 regular) under the action of Poincare return map in the 
domain where lines h = const are of the type A. rrn/M = 0.005, Ji = 1.46v^, J2 = 1.4^2. 



Fig. 12 straight line 2) gives for the lower bound- 
ary of this domain the scaling h ~ — e°-^^. The 
lower boundary of this domain was defined by the 
condition that there are no visible islands in the 
numerical experiment above this boundary. 

In the intermediate region 
—K2^hie^/^ , there are both regular and chaotic 
trajectories, and apparently the measures of sets of 
corresponding phase points are of the same order 

£2/5. 

In any annulus 7 < < /3 < with values h 
of the type A the dynamics is described by KAM 
theory (here /3, 7 are constants, i.e. values indepen- 
dent of £,7 < (3). Such an annulus is filled up to 
residual set of measure O {^/e) by invariant curves 
enclosing cylinder these curves are O {^)-c\ose 
to circles h = const. The residual set of measure 
O (a/e) is filled up to residual set of much smaller 
measure, apparently O (exp(— const/-\/£)), by in- 
variant curves that do not enclose cylinder $ 2]. 
These curves are organized into resonant " islands" . 
In this annulus, therefore, chaotic motion can oc- 



cupy only measure O (exp(— const/^/i-)). 

Let us now consider the case when value /i = is 
of the type B. Consider in cylinder $ the annulus 
a < h < Q such that all values h for this annulus 
are of type B (here a is a constant, a does not 
depend on £). In this annulus map 11 is continu- 
ous, but it is not smooth. Singularities of 11 are 
described in Sec. IIIII KAM-theory is not applica- 
ble to such a map. Trajectories of several initial 
points under the action of II in this case are shown 
in Fig. 3. Instead of invariant curves, some in- 
variant sets of complicated structure appear, that 
apparently have non-zero measure. (The structure 
of one of these sets under magnification is shown in 
Fig. 4; for similar structure for discontinuous map 
see Q and references there). Numerics show that, 
similarly to the smooth case, in this singular case 
there are several regions (see Figs. 5-6): 1) the 
region where there are no visible islands; in this re- 
gion apparently almost all trajectories are chaotic 
and correspond to escape to infinity; 2) the re- 
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Figure 2: Maximal values of h for the invariant curves enclosing cylinder $ and minimal value of h such that 
there are no visible islands at bigger values of h. The line ft = is of the type A. 



gion where there are islands, there are phase points 
corresponding to quite fast escape to infinity and 
there are also phase points outside the islands that 
do not escape for very long time (10^ iterations of 
map n do not lead to escape); about these latter 
phase points it is not clear if they are captured for- 
ever or not; 3) the region were there are island and 
all phase points outside the islands do not escape 
for very long time; again about these latter phase 
points it is not clear if they are captured forever or 
not. 

One of islands inside the domain where lines h = 
const are of the type B is shown in Fig. 7. For all 

initial points inside such an island the light shell on 
the period of its motion passes through the heavy 
shell with non-zero relative velocity and the num- 
ber of passages of the shells through each other is 
the same. Therefore, the Poincare return map on 
such an island is an analytic map and KAM-theory 
is applicable for description of dynamics inside this 
island. 



Numeric simulations give for the lower boundary 
of the domain where there are no visible islands the 
scaling h ~ —gO-ss ^ggg pjg^ 



V. DYNAMICS OF MOTION AND 
POINCARE RETURN MAP FOR THE NON 
RESTRICTED PROBLEM. 



In the non restricted problem the motion of 
shells is described by Hamiltonian H = Hi + H2. 

In this problem H is the total energy of the system, 
and it is an integral of motion. Here 



Hi 
Vi 



Vi, 



2mi 
Vio, 

Vio - Gm2mi/ri, 



ri < r-2, 
ri > 7-2, 



(7) 
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Figure 3: Trajectories of seven initial points under the action of Poincare return map in the domain where 
hnes h =const are of the type B and there are no visible escaping trajectories. mi/M — 0.0003, Ji — 1.4\/2, 
J2 = 1.2^2, Hi = -0.09. 



where 



and 



with 



G(M + TOi/2)mi Jfmi 



H2 = 
V2 = 



V20, 
V20 



V2, 

if r2 < ri, 
' Gmim2/r2, if ^2 > ri, 



Fan = 



G(M + TO2/2)to2 .7^1712 



r2 



24 



(8) 



Here pi is the momentum canonically conjugated 
to radius ri , and p2 is the momentum canonically 
conjugated to radius r2. Therefore, in the nonre- 
stricted problem motion of the shells is described 
by a Hamiltonian system with two degrees of free- 
dom. 



We have constructed Poincare return map for 
the non restricted problem. The total energy 
of system H is fixed and the shells with num- 
bers "1" and "2" are now in equal rights. We 
have constructed the Poincare return map in the 
same way as in section IIIII At the time mo- 
ments t, when the values of r2 have local minima 
on the trajectory, we mark on the cylinder <f> = 
{{h,(f)) : —00 < h < 00, (p mod 27r} the values of 
the shell "2" energy H2, and the shell "1" mean 
anomaly li : {H2,li) — {h,(j)). We present one 
such trajectory for equal masses mi = m2 = 0.007 
in FiglSI One see here chaotic diffusion and a com- 
plicated structure of island. 



VI. CONCLUSION. 

The shell approximation is often used for de- 
scription of spherical stellar clusters dynamics. We 
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Figure 4: The structure of invariant set from Fig. 3 under magnification. 



have considered dynamics of two shells in the case 
when one can neglect the influence of one (" light" ) 
shell onto the motion of another ("heavy") shell. 
It is demonstrated that dynamics of the light shell 
at small enough absolute value of negative energy 
is chaotic and leads to escape of the shell to infin- 
ity. There are two types of parameter values for 
the system. For one type the Poincare return map 
is analytic in the domain of small negative energies 
of the light shell. Therefore, KAM-theory is appli- 
cable for the description of the dynamics provided 
that the mass of the heavy shell is much smaller 
than the mass of the central body. In this case we 
have described the structure of the phase space and 
have obtained the scaling laws for the measure of 
the domain of chaotic motion and for the minimal 



energy of the light shell sufficient for escape to the 
infinity. For another type of parameter values, the 
Poincare return map has singularities in the do- 
main of small negative energies of the light shell. 
In this case KAM-theory is not applicable and it 
looks like currently there is no general theory which 
would allow to deal with this situation. Numeri- 
cal calculation show coexistence of invariant sets 
of complicated structure and chaotic trajectories 
with final escape to infinity. 
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Figure 7: Stability island inside the domain where lines h = const are of type B. 
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Figure 8: Trajectory of one initial point under the action of Poincare return map for the nonrestricted problem. 
Parameters and initial conditions are: mi = 1x12 — 0.007, Ji = 1.5^, J2 = 1.51\/2, m/mi = -0.05, H2/m2 = 
—0.05, ri = 5.1, r2 ~ 5.0, the initial direction of the shells motion is inwards. 



